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Abstract 

We perform global analysis of the solar neutrino data including the day and night 
spectra of events at SNO. In the context of two active neutrino mixing, the best fit of 
the data is provided by the LMA MSW solution with Am^ = 6.15 • 10"^ eV^, tan^ 9 = 
0.41, fs = 1.05, where fs is the boron neutrino flux in units of the corresponding 
flux in the Standard Solar Model (SSM). At 3a" level we find the following upper 
bounds: tan^ 6* < 0.84 and Am^ < 3.6 • 10"^ eV^. From Icr-interval we expect the 
day-night asymmetries of the charged current and electron scattering events to be: 
^gc^ = 3.9+1^% and = 2.1+^;^%. The only other solution which appears at 

3(T-level is the VAC solution with Am^ = 4.5 • 10"^° eV^, tan^ 9 = 2.1 and fs = 0.75. 
The best fit point in the LOW region, with Am^ = 0.93 • 10"^ eV^ and tan^ 9 = 0.64, 
is accepted at 99.95% (3.5cj) C.L. . The least point from the SMA solution region, 
with Am^ = 4.6-10~^ eV^ and tan^ 9 = 5T0~^, could be accepted at 5.5(j-level only. In 
the three neutrino context the influence of ^13 is studied. We find that with increase of 
^13 the LMA best fit point shifts to larger Am^, mixing angle is practically unchanged, 
and the quality of the fit becomes worse. The fits of LOW and SMA slightly improve. 
Predictions for KamLAND experiment (total rates, spectrum distortion) have been 
calculated. 



Pacs numbers: 14.60.Lm 14.60.Pq 95.85.Ry 26.65.+t 



1 Introduction 



The SNO data |l|, ^, ^, ^, |^ is the breakthrough in a long story of the solar neutrino problem. 
With high confidence level we can claim that solar neutrinos undergo the fiavor conversion 

z/e— >z/^,z/^ or/and z7^, z/^. (1) 

Moreover, non-electron neutrinos compose larger part of the solar neutrino fiux at high 
energies (also partial conversion to sterile neutrinos is not excluded). The main issue now 
is to identify the mechanism of neutrino conversion. 

There are several important pieces of new information from the recent SNO publication 
iii: 

1. Measurements of the energy spectra with low threshold (as well as angular distribu- 
tion) of events allow to extract information on the neutrino neutral current (NC), charged 
current (CC), as well as electron scattering (ES) event rates. In particular, in assumption 
of absence of distortion, one gets for the ratio of the NC/CC event rates: 

NC 

^ = 2.9 ±0.4 (2) 

which deviates from 1 by about 5a. 

2. Measurements of the day and night energy spectra allow to find the D-N asymmetries 
of different classes of events. Under constraint that total fiux has no D-N asymmetry one 
gets for CC event rate |Q 

Ag^ = 7.0±4.9lii%. (3) 

3. No substantial distortion of the neutrino energy spectrum has been found. 

4. Solutions of the solar neutrino problem based on pure active - sterile conversion, 
Ue ^ i^s, are strongly disfavored. 

These results further confirm earlier indications of appearance from compari- 

son of fiuxes determined from the charged current event rate in the SNO detector H , and 
the i/e— scattering event rate obtained by the Super-Kamiokande (SK) collaboration |P, |^, |^. 

Implications of the new SNO results for different solutions (see |^, ^ for earlier studies) 
can be obtained immediately by comparison of the results (0, |^) with predictions from the 
best fit points of different solutions fTT], 0, [T|, 0, 0, 0, ITTf. In particular, for the LMA 
solution the best fit prediction NC/CC = 3.3 (for lower threshold) |jl5[ is slightly higher 
than @). So, new results should move the best fit point to larger values of mixing angles. 
The expected day-night asymmetry in the best fit point, ~ 6%, is well within the interval 
d^). Clearly new data further favor this solution. 

For LOW solution: NC/CC = 2.4 [|l|] in the best fit point, which is Icr (experimental) 



lower than the central SNO value. The expected asymmetry was lower than (^, therefore 
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this solution is somewhat less favored, and SNO tends to shift the allowed region to smaller 
values of 9 which correspond to smaller survival probability. 



Implications of new SNO results have been studied in [0, |T^, ^ |22|. In this paper 
we continue this study. We perform global analysis of all available data including the SNO 
day and night energy spectra of events, and the latest data from Super-Kamiokande and 
SAGE. We identify the most plausible solutions and study their properties. 

The paper is organized as follows: In Section 2 we describe features of our analysis. 
In Section 3 we present results of the test, and construct the pull-off diagrams for 
various observables. In Section 4 we determine the regions of solutions and describe their 
properties. In Section 5 we consider the effect of 6'i3 on the solutions. In Section 6 we 
study the predictions to KamLAND for the parameters given by the found solutions. The 
conclusion is given in Section 7. 



2 Global analysis of the solar neutrino data 

In this section we describe main ingredients of our analysis. We follow the procedure of 
analysis developed in previous publications @, 



io,na,M, 



2.1 Input data 

We use the following set of the experimental results: 
1) Three rates (3 d.o.f.): 

(i) the y4r— production rate Qat measured by the Homestake experiment 

(ii) the Ge— production rate, Qce, from SAGE [p5[] , 

(iii) the combined Ge— production rate from GALLEX and GNO |2^ . 



2) The zenith-spectra measured by Super-Kamiokande during 1496 days of operation. 
The data consists of 8 energy bins with 7 zenith angle bins in each, except for the first and 
last energy bins, which makes 44 data points. We use the experimental errors given in [0 
and we treat the correlation of systematic uncertainties as in |T^. Following the procedure 



outlined in we do not include the total rate of events in the SK detector, which is not 



independent from the spectral data. 

3) From SNO, we use the day and the night energy spectra of all events [Q. We follow 
procedure described in 0. Additional information on how to treat the systematic uncer- 
tainties was given by |^ . 



Altogether there are 81 data points. 
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2.2 Neutrino Fluxes 



All solar neutrino fluxes (but the boron neutrino flux) are taken according to SSM BP2000 
We use the boron neutrino flux as free parameter. We deflne dimensionless parameter 

Fb 



28 



pSSM ■ 



5.05 • 10^ cm 



where the SSM boron neutrino flux is taken to be Fg^^^ 
/iep— neutrino flux we take flxed value Fh^p = 9.3 x 10^ cm~^ s^^ |28|, E9 



(4) 



For the 



2.3 Neutrino mixing and conversion 

We perform analysis of data in terms of mixing of two flavor neutrinos and three flavor 
neutrinos. 

In the case of two neutrinos there are two oscillation parameters: the mass squared 
difference, Am^, and the mixing parameter tan^ ^. So, we have three flt parameters: Am?, 
tan^ ^, /s, and therefore 81(data points) - 3 = 78 d.o.f. 

In the case of three neutrino mixing we adopt the mass scheme which explains the solar 
and the atmospheric neutrino data. In this scheme the mass eigenstates vi and V2 are 
splitted by the solar Am^g, whereas the third mass eigenstate, 1/3, is separated by larger 
mass split related to the atmospheric Am^3. Matter effect influences very weakly mixing 
(flavor content) of the third mass eigenstate. The effect of third neutrino is reduced then 
to the averaged vacuum oscillations. In this case, the survival probability equals 

Pee = COS^^13Pi'^+sin^^l3, (5) 

where sin ^13 = Ue^ describes the mixing of electron neutrino in the third mass eigenstate 
and Pg^g-* is the two neutrino oscillation probability characterized by tan^ 6'i2, Am^2 ^^^^1 the 
effective matter potential reduced by factor cos^ 6*13 (see e.g. |]3y, |3T| for previous studies). 

In general, in the three neutrino case the flt parameters are tan^^i2, Am^g; sin6'i3 and 
/s. However, here for illustrative purpose we take flxed value of ^13 near its upper bound. 
So, the number of degrees of freedom is the same as in the two neutrino case. 



2.4 Statistical analysis 

We perform the test of various oscillation solutions by calculating 

2 _ 2 I 2 I 2 (p\ 

X global — Xrate ' XsK ' XsNO^ I" J 

where xlatey Xsk Xsno contributions from the total rates, from the Super- 

Kamiokande zenith spectra and the SNO day and night spectra correspondingly. Each of 
the entries in Eq.(^ is the function of three parameters {Am?, tan^ ^, Jb)- 
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Some details of treatment of the systematic errors are given in the Appendix. 

The uncertainties of contributions from different components of the solar neutrino flux 
(pp-, Be-, B- etc.) to Ge-production rate due to uncertainties of the cross-section of the 
detection reaction z/g — Ga correlate. Similarly, uncertainties of contributions to Ar produc- 
tion rate due to uncertainty in v^ — Cl cross-section correlate. Following we have taken 
into account these correlations. 



2.5 Cross-checks. Comparison with other analysis 

We have checked our results performing two additional fits: 

1) To check our treatment of the SK data we have performed global analysis taking from 
SNO only the CC-rate. That corresponds to the analysis done in |^ ^. We get very good 
agreement of the results. 

2) To check our treatment of the latest SNO data we have performed analysis using the 
day and night spectra from SNO, as in We have reproduced the results of paper with 
a good accuracy. 

Our input set of the data differs from that used in other analyses: We include more 
complete and up-dated information. SNO 0] uses the SK day and night spectra measured 
after 1258 days. In contrast, we use preliminary SK zenith spectra measured during 1496 
days. In the NC/CC ratio and the D-N asymmetry at SNO where included in the 



analysis. The analysis done by Barger et al [|T8| uses the same data set we do. 



3 test 

In this section we describe the results of fit for two neutrino mixing. 

In Table |1] we show the best fit values of parameters Am^, tan^ 9, fs for different 
solutions of the solar neutrino problem. We also give the corresponding values of Xmm ^-^^ 
the goodness of the fit. 

The absolute minimum, = 65.2 for 78 d.o.f., is in the LMA region. The vacuum 
oscillation is the next best. It, however, requires ~ 30% lower boron neutrino flux. The 
LOW solution has slightly higher x^- The SMA gives a very bad fit. 

In order to check the quality of the fits we have calculated predictions for the available 
observables in the best fit points of the global solutions (see Table 1). Using these predictions 
we have constructed the "puU-ofF' diagrams (fig. |T|) which show deviations, D^, of the 
predicted values of observables K from the central experimental values expressed in the la 

K K 

Dk = — M. X = QArj Qgc, NC/CC, Ryej Adn^^dn- iV 
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Table 1: Best-fit values of the parameters Am^, tan^ 6' and fs, as well as the minimum 
and the corresponding g.o.f. for various global solutions. The number of degrees of freedom 
is 78. 



Solution Am^/eV^ tan^ 6* fs Xmin g-O-f- 

LMA 6.15 X 10^^ (L4l LOS 65^2 85%" 

VAC 4.5 X 10-1° 2.1 0.749 74.9 58% 

LOW 0.93 X 10-^ 0.64 0.908 77.6 49% 

SMA 4.6 X 10"*^ 0.5 x lO^^ 0.57 99.7 4.9% 



Here ax is the one sigma standard deviation for a given observable K. R^e is the reduced 
total rate of events at SK. We take the experimental errors only: ax = cr^^. 

According to Fig. |l| only the LMA solution does not have strong deviations of predic- 
tions from the experimental results. LOW and VAC solutions give worse fit to the data. 



4 Parameters of solutions 

We define the solution regions by constructing the contours of constant (68, 90, 95, 99, 
99.73 %) confidence level with respect of the absolute minimum in the LMA region. Follow- 
ing the same procedure as in [|1^], for each point in the Am^, tan^ 9 plane we find minimal 
value Xmin (^"^^ ! tan^ 6') varying fs- We define the contours of constant confidence level by 
the condition 

xLn(Am2, tan^ 6) = xLniLMA) + Ax' , (8) 

where x^^^{LM A) = 65.2 is the absolute minimum in the LMA region and Ax^ is taken 
for two degrees of freedom. 



4.1 LMA 

Recent SNO data further favors the LMA MSW solution (see e.g. In the best fit point 

we get 

Am' = 6.15 ■ 10^^ eV', tan^ 9 = 0.41, fs = 1.05. (9) 

The value of Am^ is slightly higher than that found in the SNO analysis and higher than 
in our previous analysis [1^. The shift is mainly due to updated SK results which show 
smaller D-N asymmetry than before. Large SNO asymmetry which would push Am^ to 
smaller values is still statistically insignificant. The mixing angle is shifted to larger values 
(in comparison with previous analysis) due to smaller ratio of the NC/CC event rates. The 
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boron neutrino flux is 5% higher than central value in the SSM: Fb = Jb ■ -Ff"^^^ = 5.32 • 10^ 
cm~^ being however within la deviation and well in agreement with SNO measurements. 

The CL contours (see flg. ^ shrink substantially as compared with previous determina- 
tion fn], ll, m 0. 



From flg. ^ we flnd the following bounds on oscillations parameters: 

1) Am^ is rather sharply restricted from below by the day-night asymmetry at SK: 
Am2 > 2.3 ■ 10-5 eV2 at 99.73% CL. . 

2) The upper limits on Am^ for different confldence levels equal: 

1.2 X 10-^ eV^ 68.27% CL. 
Am^ <{ 1.9 X 10-^ eV^ 95% CL. . (10) 
3.6 X 10-^ eV^ 99.73% CL. 




All these limits are stronger than the CHOOZ ||3^ bound which appears for maximal mixing 
at Am2 ~ 8 X 10-^eV^ 

3) The upper limit on mixing angle becomes substantially stronger than before: 

68.27% CL. 

tan^9< { 0.65 95% CL. (11) 
99.73% CL. 

Maximal mixing is allowed at ~ 4cr level for Am^ = (5 — 7) ■ 10"^ eV^. 

Notice that the SNO data alone exclude maximal mixing at about 3a: the data deter- 
mine now rather precisely the NC/CC ratio which is directly related to sin^ 9. Also observed 
Germanium production rate as well as Argon production rate disfavor maximal mixing (see 
flg. I and I). 

So, now we have strong evidence that solar neutrino mixing signiflcantly deviates from 
maximal value. One can introduce the deviation parameter |^| 



e = l-2sm^9. (12) 

From our analysis we get 

e > 0.08, (3a). (13) 

That is, at 3a: e > sin^ 9^, where 9^ is the Cabibbo angle. This result has important 
theoretical implications. 

4) lower limit on mixing : 

tan2^> 0.23, 99% CL.. (14) 

is changed weakly. 
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In fig. we show the grids of predicted values for various observables. 

According to the pull-off diagram and figs. ^||, the LMA solution reproduces observables 
at ~ 1(7 or better. The largest deviation is for the Ar— production rate: the solution predicts 
1.6(7 larger rate than the Homestake result. 

The best fit point value and 3a interval for Ge production rate equal 

Qce = 70.5 SNU, Qce = (63 - 84:)SNU, 3a. (15) 

Notice that at maximal mixing Qce < 63 SNU which is 2a away from the combined exper- 
imental result. 



4.2 VAC 

In the best fit point we get = 74.9 and: 

x\VAC)-x\LMA) = 9.7. (16) 

So, this solution is accepted at 3a level. Notice that the solution appears in the dark side 
of the parameter space which means that some matter effect is present. This solution was 
"discovered" in 1998 and its properties have already been described in the literature. Clearly 
it does not predict any day-night asymmetry. The solution requires rather low (1.6(7) Boron 
neutrino flux and gives rather poor description of rates (see flg. |1|). In particular, 2.7a 
higher Ar-production rate and 2.6a lower NC/CC ratio are predicted. Imposing the SSM 
restriction on this flux leads to exclusion of this VAC solution at 3a level. 



4.3 Any chance for SMA? 

We flnd that the best flt point from the SMA region has = 99.8. For the difference of 
we have: 

X^iSMA) - x^iLMA) = 34.5. (17) 

That is, SMA is accepted at 5.5cr only. Moreover, the solution requires about 3a lower boron 
neutrino flux than in the SSM. It predicts negative Day-Night asymmetry: = —0.93%. 

Our results are in qualitative agreement with those in ||18|, where even larger Ax^ has 
been obtained. 

We find that the x^ increases weakly with tan^ 6 up to tan^ 6 = 1.5-10^^, where x^ ~ 105. 

Is SMA excluded? We find that very bad fit is due to latest SNO measurements of 
day and night spectra. We have checked that the analysis of the same set of data but CC 
rate from SNO only (2001 year) instead of spectrum leads to the best fit values Am^ = 
4.8 ■ 10-6 eV^ tan^^ = 3.9 • 10"^ and = and x^iSMA) - x^{LMA) = 11 in a good 
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agreement with results of similar analysis in [§] . Since CC SNO data are in a good agreement 
with new NC/CC result, just using the NC/CC does not produce substantial change of 
quality of the SMA fit |0. So it is the spectral data which give large contribution to x^- 

The SMA solution with very small mixing provides rather good description of the SK 
data: the rate and spectra. The (reduced) rate R = OBS / SSM of the ES events can be 
written as 

[ES] = fB[Pee{l-r) + r], (18) 

where r 0.16 is the ratio of z/^ — e to z/g ~ e cross-sections. Taking Res = 0.45 and 
fs = 0.58 we find the effective survival probability: Pge = 0.73. Then, for reduced CC event 
rate we get [CC] = fsPee = 0.425 - close to the ES rate, and moreover, 

NC/CC ^ 1/Pee = 1.37 (19) 

which is substantially smaller then the observed quantity (H). So, one predicts in this case 
a suppressed contribution of the NC events to the total rates. Correspondingly, significant 
distortion of the energy spectrum of events is expected with smaller than observed rate at 
low energies and higher rate at high energies. 

This problem with SNO could be avoided for larger mixing: tan^ 6 > 1.5 • 10^'^ (in fact, 
imposing the SSM restriction on the boron neutrino flux leads to the shift of the best fit 
point to larger 6). In this case, however serious problems with SK data appear, namely, 
with spectrum distortion and zenith angle distribution. Strong day-night asymmetry is 
predicted for the Earth core-crossing bin. Previous analysis which used SK day and night 
spectra could not realize the latter problem. 

Notice that the SNO data alone do not disfavor SMA with large tan^ 9 = (1.5-2) ■ 10"^. 
This region, however is strongly disfavored by SK. 

Zenith angle distribution can give a decisive check of the SMA solution. The SNO night 
data could be divided into two bins: "mantle" and "core". Concentration of the night 



excess of rate in the core bin ^J] due to parametric enhancement of oscillations for the core 
crossing trajectories would be the evidence of the SMA solution with relatively large 
mixing: tan^^ = (1.5 — 2) • 10~^. However, the SK zenith spectra do not show any excess 
of the "core" bin rate which testify against this possibility. 

Probably some unknown systematics could improve the SMA fit. Otherwise, this solu- 
tion is practically excluded. 

4.4 LOW starts to disappear? 

In the best fit point we get = 78.9, so that 

x\LOW) - x\LMA) = 12.4 (20) 
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which is shghtly beyond the 3a range. In contrast with other analyses, LOW does not appear 
at 3a leveL Notice that in the SNO analysis |Q the LOW solution exists marginally at 3a 
level. Inclusion of the SK data which contain information about zenith angle distribution 
(zenith spectra) worsen the fit (this effect has also been observed in |I3| ). 



The LOW solution gives rather poor fit of total rates. In the best fit point we get 2.1a 
larger Ar— production rate and 1.2a lower Ge— production rate. For the day-night asym- 
metry of the CC-events we predict = 3.5% and for ES events: A'f)% = 2.7%. 



5 Three neutrino mixing: effect of ^13 

Results of the global analysis in the three neutrino context are shown in fig. |^. To illustrate 
the effect of third neutrino we use the three neutrino survival probability (^ for fixed value 
sin^ ^13 = 0.04 near the upper bound from the CHOOZ experiment The number of 

degrees of freedom is the same as in the previous analysis. 

We find the best fit point: 

Am?2 = 6.7 • 10" tan^ 9 = 0.41, /b = 1.09 (21) 

with = 66.2. The best fit value of Am^g is slightly higher than that in the two neutrino 
case, whereas the mixing angle is unchanged. The solution requires slightly higher value of 
the boron neutrino fiux. The changes are rather small, however, as a tendency, we find that 
with increase of Oi^ the fit becomes worser in comparison with 2i/— case. For sin^ = 0.04 
we get Ax^ = 1-0. 

In fig. 1^ we show the contours of constant confidence level constructed with respect to 
the best fit point (pi]). The contours changed weakly for low mass values Am^g < 10~^ eV^ 
and there are significant changes for Am^g > 10~^ eV^. In particular, the 3a upper bound 
on Amfg is Am^2 < 5.8 • 10~^ eV^; the lower 3a bound on mixing: tan^ Ou = 0.18 (compare 
with numbers in the Table 1). Notice, however, that changes are substantially weaker if the 
contours are constructed with respect to the absolute minimum for ^13 = (^. 

The changes can be easily understood from the following analytical consideration. 

The contribution of the last term in the probability (|^) is negligible: for largest possible 
value of 6*13 it is below 0.5%. So, we can safely use approximation: 

Pee ~ COS^ ^IS^i? ~ (1 - 2 siu^ ^^13)Pie^- (22) 

Mainly the effect of ^13 is reduced to overall suppression of the survival probability. The 
suppression factor can be as small as 0.90 - 0.92. 

In the fit with the free boron neutrino fiux, the observables at high energies ( > 5 MeV) 
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are determined by the following reduced rates 



cc 

[CC] = ^^55M= fB cos' e^sP!:!\ 

IES^^E§^^= fB[cos'0,sPiHl-r)+r]. (23) 

As far as the fit of experimental data on CC-events are concerned (SNO, SK, and partly, 
Homestake), the effects of ^13 is simply reduced to renormalization of the boron neutrino 
flux: 

fs - (24) 
cos* 

without change of the oscillation parameters Ami2 and 612. The dependence of the pa- 
rameters on 6*13 appears via the ratios of rates, which do not depend on fs- From (^31) we 
find 



ee 



cos'' 9 



13 



[ES] 



p(2)- 
lee 



(26) 



[CC] ' 

So, the effect of 6*13 is equivalent to decrease of the ratios [NC]/[CC] and [ES]/[CC]. Accord- 
ing to fig. 1^, this shifts the allowed regions to larger Ami2 and 6*12. 

For low energy measurements (Gallium experiments), sensitive to the pp-neutrino flux, 
which is known rather well, the increase of ^13 should be compensated by increase of the 
survival probability. This may occur due to increase of Ami2 or /and decrease of tan^^i2. 

For Am^2 < 10"'* eV^ the boron neutrino spectrum is in the bottom of suppression pit 
and the low energy neutrinos are on the adiabatic edge. In the fit, the increase of ^13 is 
compensated by the increase of fs and Am^2- For Am^g > 10~^ eV^, the spectrum is in the 
region where conversion is determined mainly by averaged vacuum oscillations with some 
matter corrections: Pjg-* ~ (1 — 0.5 sin^ 2^12). The dependence on Amf2 is very weak which 
explains substantial enlargement of the allowed region to large values of Am^2- The effect 
of 6'i3 can be compensated by decrease of 612 which explains expansion of the region toward 
smaller tan^ 612. 

For LOW solution increase of ^13 leads to improvement of the fit, so that this solution 
appears (for sin^ ^13 = 0.4) at 3cr level with respect to best fit point (0). Also for the SMA 
solution the effect of ^13 leads to slight improvement of the fit. 
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6 Predictions for KamLAND 



Next step in developments will be probably related to operation of the KamLAND exper- 
iment Both total rate of events above the effective threshold T^ff and the energy 



spectrum of events will be measured. 

We characterize the effect of the oscillation disappearance by the ratio of the total 
number of events with visible energy above T^ff: 

RKam = ^[ I [ dEdT'dTj2F,P,^f{T,T'), (27) 

Nq JT^ff JT' JE dT 

where Fi is the flux from i reactor, Pj is the survival probability for neutrinos from i reactor, 
a is the cross-section of the detection reaction, f{T,T') is the energy resolution. Nq is the 
rate without oscillations {Pi = 1). 



In our calculations we used the energy spectra of reactor neutrinos from p9|, E0[. The 



differential cross-section of the p + V^ ^ n + e'^ reaction, is taken from ^J]. The parameters 
of the 16 nuclear reactors, maximal thermal power, distance from the reactor to the detector, 
etc., are given in We used the Gaussian form for the energy resolution function /(T, T') 



with a/E = 5%/ ^ E{MeV), and Tg// = 2.6 MeV as the threshold for the visible energy 
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In fig. H we show the contours of constant suppression factor in the Am^ — tan^ 6 plot. 
In the best fit point 

RKam = 0.65 , (28) 

and in the la region: Rxam = 0.4 — 0.7. 

Notice that the best fit point is in the range of lowest sensitivity of the total rate on 
tan^6'. If e.g. RKam is measured with 8% accuracy which would correspond to Rxam = 
0.65 ± 0.05, we get from the fig. || that any mixing in the interval tan^ ^ = 0.12 — 1.0 is 
allowed. 

The suppression factor strongly depends on Am^ in the range below the best fit point 
and this dependence is very weak for Am^ > 10"^ eV^. No bound on Am^ from the allowed 
region can be obtained by measurements of the total rate. 

The distortion of the visible energy spectrum depends on Am^ strongly. In fig. ^ we 
show the spectrum for different values of Am^. There is a shift of maximum to large E with 
increase of Am^. For the best fit value of Am^ the maximum is at E' ~ 3.5 MeV. The most 
profound effect of oscillations is the suppression of rate at high energies. For instance, for 
E ^ 5 MeV the suppression factor is smaller than 1/2. 
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7 Conclusions 



We find that the LMA MSW solution with parameters Am^ ~ 6.15 • 10~^ eV^ and tan^ 9 = 
0.41 gives the best fit to the data. The solution reproduces well the zenith spectrum mea- 
sured by SK and the day and night spectra at SNO. It is in a very good agreement with 
SSM fiux of the boron neutrino: Jb = 1.05. 

The recent SNO results together with zenith spectra results from SK slightly shifted 
the best fit point to larger Am^ and 6. At the same time the allowed regions of oscillation 
parameters shrunk, leading to important, and statistically significant, upper bounds on 
mixing angle and Am^. Now we have strong evidence that "solar" mixing is non-maximal, 
and moreover, deviation from maximal mixing is rather large. We find that QuasiVacuum 
oscillation solution with Am^ = 4.5 ■ 10~^° eV^ and mixing in the dark side is the only other 
solution accepted at 3a level, provided that the boron neutrino fiux is about 30% below the 
SSM value. 

The LOW solution is accepted at slightly higher than 3a level and it reappears at 3a 
level if 6^13 is included. 

The SMA solution gives very bad fit of the data especially the SNO spectra predicting 
rather small contribution of the NC events in comparison with CC events. 

We find that 613 produces rather small effect on the solutions even with new high 
statistics data. As a tendency we see that inclusion of the ^13 effect worses fit of the data 
in the LMA region, and shifts the best fit point to larger Am^2- 

We have found predictions for the KamLAND experiment: in the best fit point one 
expects the suppression factor for total signal ~ 0.6 — 0.7 and the spectrum distortion with 
substantial suppression in the high energy part. 
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Appendix 



The systematics uncertainties were treated according to p^. Writing the total counting 



rate in SNO as a sum over different contributions and different spectral bins, we have: 
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Rj = ^ Rij , (29) 



i=l,5 



where the index j stands for the different spectral bins and i runs over the five contributions 
to the SNO data (CC, NC, ES, neutron background and low energy background). We assume 
that all systematic uncertainties of the SNO result are the uncertainties in the theoretical 
prediction. These uncertainties can be written in terms of the systematics uncertainties of 
the input parameters of experiment (Xk): 

We take the systematic uncertainties from The different systematic uncertainties 
are added in quadrature. 



Eq.(^) can be written in terms of the different contributions Rij ( p9D as: 

(A/nXfc)^ (31) 

ii=l,5 42=1,5 fc=l,14 

where we have introduced the parameters 

dlnRi i , , 

These parameters are numerically estimated by changing the response function of the 
detector through changes in the parameters X^. 
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Figure 1: Pull-off diagrams for global solutions. Shown are deviations of predictions from 
experimentally measured values for the Ar— production rate, Ge— production rate, SK rate, 
the day-night asymmetries at SK and SNO. The puU-offs are expressed in the units of 1 
standard deviation, la. 
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Figure 2: The global LMA MSW solution. The boron neutrino flux is considered as free 
parameter. The best fit point is marked by a star. The allowed regions are shown at la, 
90% C.L., 2a, 99% C.L. and 3a. 
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Figure 3: Lines of constant Ge-production rate (number at the curves in SNU) in the LMA 
region. In the best fit point: Rce = 70.5 SNU. 
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Figure 4: The same as in Fig ^, but for the Ar-production rate. In the best fit point: 
Rat = 2.95 SNU. The dependence of fs on oscillation parameters is taken into account. 
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Figure 5: Lines of constant NC/CC ratio in the LMA region. In the best fit point: 
NC/CC = 3.15. 
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Figure 7: Global LMA solution for sin^ ^13 = 0.04. The boron neutrino flux is considered 
as free parameter. The best fit point is marked by a star. The allowed regions are shown 
at la, 90% C.L., 2a, 99% C.L. and 3a. 
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Figure 8: Lines of constant total suppression at KamLAND. In the best fit point: RKam = 
0.65. 
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Figure 9: Spectral distortion for three different values on Am^, including the best fit point 
of our analysis. 
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